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ABSTRACT 

We present quasi-simultaneous, multi-epoch radio and X-ray measurements of Holm¬ 
berg II X-1 using the European VLBI Network (EVN), the Karl G. Jansky Very 
Large Array (VLA), and the Chandra and Swift X-ray telescopes. The X-ray data 
show apparently hard spectra with steady X-ray luminosities 4 months apart from 
each other. In the high-resolution EVN radio observations, we have detected an ex¬ 
tended milli-arcsecond scale source with unboosted radio emission. The source emits 
non-thermal, likely optically thin synchrotron emission and its morphology is consis¬ 
tent with a jet ejection. The 9-GHz VLA data show an arcsecond-scale triple structure 
of Holmberg H X-1 similar to that seen at lower frequencies. However, we find that the 
central ejection has faded by at least a factor of 7.3 over 1.5 years. We estimate the 
dynamical age of the ejection to be higher than 2.1 years. We show that such a rapid 
cooling can be explained with simple adiabatic expansion losses. These properties of 
Holmberg H X-1 imply that ULX radio bubbles may be inflated by ejecta instead of 
self-absorbed compact jets. 

Key words: accretion, accretion discs - black hole physics - X-rays: binaries 


1 INTRODUCTION 


Ultraluminous X-ray sources (ULXs) host accreting compact 
objects, whose luminosity {Lx > 3 x 10^® erg/s) exceeds 
the Eddington limit of a 2 O-M 0 black hole (BH). ULXs are 
a mixed bag o f sources. Some may host super-Eddington 


neutron stars (|Bachetti_etaL 2014 ) or stellar-ma ss {M = 
3—20 M 0 ) BHs ( Liu et al.ll2013l i^ Motch et al.ll2014l ). Some of 


BHs ( 

Farrell et al.ll2009l: 

Fens & Kaaretll2010l:lPasham et al.l 

l20l4 

Mezcua et al.ll2015 

) and may also host recoiling super- 

massive BHs (Jonker et al.||2010|). 


Given the tentative observational evidence for the anti- 
corr elation between BH mass and host galaxy metallic- 
ity dCrowther et al.l 120101 ) and that some ULX progeni- 
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tor stars could potentially have been m assive enough (e.g. 
iGrise et'al]|20l4 iPoutanG et al.ll2013l ) to collapse directly 
to a BH, ULXs m i ght host massive stel l ar-mass (20-100 M 0 ) 
BHs llFrver[ 19 ^ Jms-pcUI et al]l2009l : ZampieH fc Roberta 
2OO9I: iBelczvnski et al.l 2010l : Mapelli fc Bress^ 20131 : 


Bachetti et al.l 20131 ). In fact, Holmberg H X-1 lies in a rel¬ 


ative l y metal-poor dwar f irregular galaxy dPrestwich et al.l 
I2OI3I : Brorbv et aI.ll2oG ) . Such a mass range has been shown 
to fit its observed X-r ay and radio properties dGoad et al.l 
|2QQ^ Cseh et ^ 2014h and may also be valid for IC 342 X-1 
dMarlowe et al.ir2014l ). 


To date, relatively few ULX radio counterparts have 
been studied, and these few show a relatively diverse na¬ 
ture. The first radio bubble, associated with NGC 5408 X - 
1, was discovered about a decade ago l|Kaaret et ahllioo^ 
although one other, associated with NGC 6946 X-1, was 
detected earlier but considered as a super- or hypernova 
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remnant dvan Dvk et alJ 1 19941) . FollowinR the discoveries 
of more radio bubbles (e.g. ICseh et al.l l2012l . and ref¬ 
erences therein), transient jets were also caught in two 
non-persistent ULXs, which in both cases were associ- 
ated with X-ray outbursts of the centra l compact object 
(|Webb et al.ll201^ : iMiddleton et al.ll201,'j l. Furthermore, in 
a few other cases, apparently steady radio emission has 
been detected on milli-arcs e cond scales (iMezcua fc Lobanov! 
I2OIII : iMezcua et alJ |2013| . l2015ll and interpreted as self- 
absorbed radio cores, which are typically seen during ra- 
diatively ine fficient accretion st ates of Galactic BH X-ray 
binaries fe.g. iFender et al.]l2009l) . 

Focusing on the large-scale environments of ULXs, jets 
have been argued to be a plausible supply mechanism re¬ 
sponsible for inflating the observe d radio bubbles (and 
shock-ionized optical bubbles) le.g. ICseh et al.ll2012l ). On 
th e other hand, dir ect evidence of jets in Holmberg II X- 
1 dCseh et alJ l2014l) has been found, despite an environ- 
ment that appears to be dominated by photoionisation 
llLehmann et alJ l2005l : lEgorov et al.11201^ . This may indi¬ 
cate that shock ionization is not a neces sary consequence 
of jet s . Shock-ionized optica l bubbles fe.g. lPakull fc Mirionil 
I2OO2I : IPakull fc Grisd l2008l ) could equally be inflated by 
powerful winds. These winds are thought to originate from 
a supercritical disk that may also give rise to geometri¬ 
cally collimated X-rays llPoutanen et al.l I200i1) . although 
calorimetry of photoionized optical bubb les argues against 
stron g ly beamed X-ray emi ssion (e.g. IPakull fc Mirionil 
I2OO2I : likaaret fc Gorb J l2009l L Moreover, direct observa- 
tional evidence for super-Eddington winds in ULX s is scarce 
ll Walton et al.l [20121 . 12OI3I : IMiddleton et al.ll2014h . Interest¬ 
ingly, a similar group of objects, such as S26 in NGG 7793 
and MQl in M83, which do not appear to have luminous 
central X-ray so urces, can nonethe l ess have more ene rgetic 
optical bubbles llPakull et al.l l20ld : ISoria et al. than 

classical ULXs. These two peculiar examples both appear 
to have radio jets (or at least hot spots) and relatively large 
cocoons. Such morphology resembles SS433, although the 
total energy content of the W50 bubble surrounding SS433 
is two orders of magnitude lower ilPakull et al.ll2Mg ). 

In the Galactic BH X-ray binaries, there is an obser¬ 
vational anticorrelation between the presence of disk winds 
and jets. Winds are observed in soft or disk-dominated ob¬ 
jects and jets are typicall y seen (at lower Eddington rates) 
in hard-state objects (e.g. Ponti et al ]|2 m 3). Broadly speak¬ 
ing, in active galactic nuclei (AGNs), near-Eddington lu¬ 
minosities are also coupled with powerful winds, whereas 
relativistic jets tend t o occur at lower accretion rates (e.g. 
iKing fc PoundsI l2015l L The presence of a collimated jet 
structure in Holmberg II X-1 coupled with a high X-ray 
luminosity might be somewhat extraordinary if one as¬ 
sumes persistent near- or super-Eddington accretion. How¬ 
ever, little is known about switching between sub- and super- 
Eddington accretion. Possibly the best examples for such 
transitions are the Galactic systems GRO J1655-40 (e.g. 
iMeieillig^ and GRS1915-I-105 (e.g. lArai et ahllioogl L Ad¬ 
ditionally, the higher BH masses inferred for sources like 
Holmberg H X-1 may imply longer characteristic timescales 
llMcHardv et al.ll2006l L 

In summary, Holmberg H X-1 is a bona fide ULX 
with a typical X -ray luminosity of ~ 10'*° erg s“* (e.g. 
iGrise et al.ll201(]| ). It is hosted by a nearby, relatively metal- 


Table 1. Summary of the observations 


Inst. 

Mode 

Frequency or 

On-source 

Observation 


or config. 

energy band 

time 

date 

EVN 

- 

1.6 GHz 

12.5 h 

14 Jan 14 

Swift 

PC 

0.3 - 10 keV 

20.5 ks 

17-22 Jan 14 

EVN 

- 

5 GHz 

12.5 h 

28 Apr 14 

Chandra 

1/8 

0.3 - 10 keV 

11.8 ks 

25 May 14 

VLA 

A 

8-10 GHz 

171 min 

25 May 14 


poor (e.g. IPrestwich et ^ 20lj) dwarf irregular galaxy at 
a distance of 3.39 Mpc ( Karache ntsev et alll2002l~l. It is as¬ 
sociated wit h an optical (e.g. Kaaret et al.l 2004h and a ra¬ 
dio bubble dMiller et al.ll200^ . The optical environment is 
photoionized wi th an average isotr opic X-ray luminosity of 
4 X 10°° erg s~* (lKaaret et al.ll2004l L The radio bubble is in¬ 
flated due to jet activity a nd the deduced av erage jet power 
of Qj ~ 2 X 10°° erg s“° ijGseh et al.ll201'3 l is comparable 
to the average X-ray luminosity. Its X-ray spectrum shows 
typical ULX feature s, like a soft excess and a spectral break 
llKaiava et al.l I2OI2I) . Based on the energetics, the X-ray 
spectral properties, and the en vironment, its mass is likely to 
be ill the range of 25-100 M© llGoad et al.ll2006l : iGseh et al.l 
I 2 OI 4 II with an accretion rate near or above the Eddington- 
limit. The high accretion rate may also be supported by 
that the jet morphology indicates a recurrent activity (see 
below). Also, the above mass range might be supported by 
the longer act ivity of the sour ce as compared to Galactic 
microquasars dCseh et al.ll2014h . 

How these massive stellar-mass BHs evolve remains 
unclear. Early BH growth is often hypothesised to pro- 
ceed via (intermittent) super-Eddington acc retion (e.g. 
IVolonteri fc Silkll2014l : IPacucci fc Eerraral 120151) . To inves¬ 
tigate if some ULXs and their observable properties can act 
as a proxy for super-Eddington BH growth in the early Uni¬ 
verse, we present quasi-simultaneous and multi-epoch obser¬ 
vations of the ULX Holmberg H X-1 (hereafter Ho II X-1), 
which are summarised in Table [T] Our primary aim is to 
study the nature an d evolution of the previously-detected 
central radio source dCseh et al.ll20lj ). and to investigate 
the link between the radio emission and the X-ray behaviour. 
In Sec. [ 2 ] we describe our observational results. In Sec.[3l we 
first discuss the X-ray spectra and the high-resolution radio 
properties of Ho H X-1, followed by our interpretation of the 
temporal evolution of the radio ejecta in Ho II X-1. Finally, 
we discuss the physical properties of the jet ejecta. 


2 OBSERVATIONS, ANALYSIS, AND 
RESULTS 

2.1 Joint Sw/t/XRT and EVN observations 

2.1.1 Swift/XRT observations 

The Swift X-ray telescope (XRT) observed Ho H X-1 (PI: 
Grise) in photon counting (PC) mode on 17, 20, and 22 Jan 
2014 starting from 13:36:59 UT, 00:48:59 UT, and 02:27:58 
UT, for useful exposures of 5.71 ks, 5.88 ks, and 8.96 ks, 
respectively. We retrieved level 2 event files that were sub- 
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jected to standard screening method^ The source spec¬ 
trum was extracted using an aperture of 20 pixels, which 
corresponds to 90% of the PSF at 1.5 keV. The back¬ 
ground was extracted using an aperture of radius 40 pix¬ 
els, centred on a J2000 location of RA=08:19:51.437 and 
Dec=-|-70:47:46.29. The auxiliary response file (ARF) was 
created with xrtmkarf, which accounts for PSF correction. 

After rebinning the spectrum to have a minimum 
of 30 counts per bin, we fit the X-ray spectrum using 
XSPEC V12.8.2 Ia rnau dl ll996l ~l and a response matrix of 
swxpcOtol2s6_2O13OlOlvO14.rm0. We fit the spectrum with 
an absorbed power law model (tbabs*tbabs*powerlaw) 
and froze the first hydrogen absorption column den- 
sity to the Galacti c value of Nh = 3.4 x 10^° cm“^ 
dPickev fc Lockm^l 19901 1. We found a good fit of x^/DoF= 
0.99 for 80 degrees of freedom (DoF), a photon index of 
F = 1.91 ± 0.08, an equivalent hydrogen absorption column 
density of Nh = (5 ± 2) x 10^° cm~^, and an unabsorbed 
flux of (6.6±0.4) X 10“^^ erg cm“^ s“^ in the 0.3-10 keV en¬ 
ergy range. The fitted power-law photon index of F = 1.91 
is consistent with a hard X-ray spectrum. The flux corre¬ 
sponds to a luminosity of (9 .1 ± 0.6) x 10^^ erg s~^ a t the 
source distance of 3.39 Mpc dKarachentsev et al.ll2002 ). We 
fit other models to the X-ray data. First a disk blackbody 
model, that gave a poor fit with a reduced > 3- A. cut¬ 
off power law model gave an adequate fit, but with a cutoff 
energy at 500 keV, indicating that a simple power law fits 
well. A combination of a disk blackbody with a power law 
resulted a 0 keV inner disk temperature, indicating that a 
simple power law fits well. _ 

Motivated by the best-fit model of lWalton et al.l ll2015ll , 
we also fit the Swift spectrum with two disc blackbody 
models (TBABS*TBABS*(DlSKBB-t-DlSKBB)). This fit resulted 
in a reduced x^/DoF=0.94, that may indicate an overcon¬ 
strained model. The first hydrogen absorption column den¬ 
sity was frozen to the Galactic value. The fit did not require 
any additional Nh on top of the Galactic one. The fit re¬ 
sulted in inner disc temperatures of Tin,i = 0.32 ±0.05 keV 
and Tin ,2 = l- 7 lo ’2 keV. The hotter component a grees 
within the errors with the result of lWalton et al.l (l2015ll and 
the cooler component is inconsi stent within the erro rs; it is 
somewhat hotter than found by I Walton et al.l ll2015ll . 


2.1.8 VLBI observations 

We carried out very long baseline interferometry (VLBI) 
observations of the central component of Ho II X-1 with 
the electronic European VLBI Network (e-EVN) under pro¬ 
gram codes EG045 and RCOOl (PI: Gseh), at 1.6 GHz 
and 5 GHz, respectively. The 1.6-GHz observation started 
on 14 Jan 2014 15:01:10 UT and had a total duration of 
about 18 h. The 5-GHz observation started on 28 Apr 2014 
13:01:10 UT and lasted for about 18 h. At 1.6 GHz, the 
participating stations were Effelsberg (Ef; Germany), Jo- 
drell Bank Lovell Telescope (Jb; United Kingdom), Medic- 
ina (Me; Italy), Noto (Nt; Italy), Onsala (On; Sweden), She- 
shan (Sh; China), Toruh (Tr; Poland), and Westerbork (Wb; 


^ http://heasarc.nasa.gov/docs/swift/analysis/ 

^ http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/- 
data/swift /xrt/index.html 
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Figure 1. The 1.6-GHz e-EVN image of Ho II X-1 that was made 
using natural weighting. The lowest contours start at three times 
the rms noise level, at ±15 /rjy beam~^. Contours are increased 
by a factor of ^/2. The peak brightness is 53 /rjy beam“^. The 
Gaussian restoring beam size is 12 mas X 9.6 mas at a major axis 
position angle of 32°. 


The Netherlands). For the 5-GHz observation, the stations 
were Ef, Jb Mark H, Nt, On, Tr, Wb, Sh, and Yebes (Ys; 
Spain). The data recording rate was 1 Gbit s“^ per station 
(except Me where it was 512 Mbit s“^) and no data were 
received from Ys. There were eight subbands, each with 16 
MHz bandwidth in both left and right circular polarisation. 

The observations were carried out in standard phase- 
referencing mode, using J0841±7053 as the phase- 
calibrator. The observations cycled continuously between 
the target and the reference source, spending 3.5 min on 
the target in every 5- min block. Th e data were calibrated in 
AIPS v3 1Decl2 (e.g. GreisenllioO^ following standard pro¬ 
cedures toiamon dlll99,5ll and imaging was performed using 
Difmap (IShenherd et al.lll994l L Antenna-based gain correc¬ 
tion factors were obtained from the phase calibrator after a 
global amplitude self-calibration step, which were fed back 
to AIPS and applied to the data. These corrections were 
typically less than 10% at 1.6 GHz and less than 20% at 5 
GHz. Fringe-fitting was then repeated taking into account 
the clean component model of the calibrator. 

At 5 GHz, we reached an rms noise level of 12 pjy 
beam“^, the Gaussian restoring beam size was 2.7 mas x2.1 
mas, and we did not detect Ho H X-1. We set a 3 -ct upper 
limit of 36 /rJy on the flux density of the source. 

At 1.6 GHz, the rms noise level was 0^.5 ^Jy beam”*^ 
and the Gaussian restoring beam size was 12 mas x 9.6 mas 
in PA 32°. Ho H X-1 was detected at a signal-to-noise ratio 
of > 10 (Fig. [ 1 ]), with a peak brightness of 53 /rjy at a posi¬ 
tion of (J2000) RA=08:19:28.9835 and Dec=±70:42:18.9966 
(accurate to within 1 mas). To check the compactness of the 
source we fit the brightness distribution with a Gaussian 
in both the image and the visibility plane. Using the AIPS 
task JMFIT for the image-plane fit, we found a fitted peak 
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brightness of 49 ± 6 pJy beam“^ and an integrated flux 
density of 84 ± 14 /rJy. The integrated flux density is con¬ 
siderably higher than the peak brightness, likely indicating 
a resolved source. 

A fit in the visibility plane using a point-source model in 
Difmap resulted in a total flux densitjQ of Si, = 65 ± 6 /rJy. 
If we fit the source with a circular Gaussian model, we find 
a total flux density of S,/ = 88 ± 11 /rjy and a full width at 
half maximum (FWHM) size of 16.4 ± 1.7 mas. A compar¬ 
ison of these flux densities also indicates that the source is 
likely resolved. The obtained size of 16.4 ± 1.7 mas exceeds 
the restoring beam size of 12 x 9.6 mas. Furth ermore, we 
estimate a resolution limit (iKovalev et ahlBoOSl l of ~3 mas 
that is well sampled by the image cell size of 1 mas. 

To test at what angular resolution the source starts 
to be resolved, we produced an image with a restricted 
(u, v) range of 0-4 MA. The resulting image had a beam 
size of 33 x 30 mas. The source had a peak brightness 
of 67 /iJy beam“^, and a fit to a circular Gaussian model 
gave a total flux density of 95±10 fj,Jy with a FWHM of 
20 ± 1.6 mas. There is an apparent increase in all of the 
fit parameters, although they are consistent with the above 
values within the errors. We conclude that the source is re¬ 
solved on VLBI scales. We set a lower limit to its angular 
size of 16.4T1.7 mas that corresponds to a projected size of 
0.26 ± 0.03 pc at a distance of Ho H X-1. A short summary 
is shown in Table [5] Should the true size be larger than 
16.4 mas, we estimate an upper limit of ~ 60 mas from our 
highest resolution VLA images. 

2.2 Joint Chandra and VLA observations 

2.2.1 Chandra observations 

We observed Ho H X-1 with the Chandra X-ray Observatory 
under ObsID 15771 (PI: Jonker) using the Advanced CGD 
Imaging Spectroscopic array (AGIS) S3 chip in FAINT mode 
and with 1/8 subarray. The observation was carried out on 
25 May 2014, started at 03:21:52 UT, and had a useful ex¬ 
posure of 11.85 ks. There were no strong background flares. 
The data were processed in a standard way and we used the 
CIAO version 4.6 and the CALDB version 4.6.1 to further 
process the level 2 event files. 

First, we checked for background flares following stan¬ 
dard procedure^. We removed strong sources and then 
checked the extracted background light curve and we found 
no strong flares during the observation. However, Ho H X- 
1 is a relatively bright X-ray source that caused a read¬ 
out streak, which we therefore removed as described in the 
Chandra science thread page^. Using wavdetect, we de¬ 
tected Ho H X-1 with a total of 6406 photons, and extracted 
its spectrum in the 0.3-10 keV range using specextract, 
with a circular aperture of radius 3 arcsec. 

Given that we employed 1 /8 subarray, providing a frame 
time of 0.4 sec, the expectecfl pileup fraction of Ho H X-1 is 
~5%. To fit the spectrum we used the Sherpa package and its 

® The error on this and the subsequent flux densities includes a 
5% systematic uncertainty. 

http://cxc.harvard.edu/ciao/threads/fllter_ltcrv/ 

® http://cxc.harvard.edu/ciao/threads/acisreadcorr/ 

® http://cxc.harvard.edu/toolkit/pimms.jsp 


built-in pile-up model, following the Sherpa thread pagefl 
First, we grouped our spectrum to have 30 counts per bin 
and then subtracted the background. We used an absorbed 
power-law model (xstbabs*xstbabs*xspowerlaw), set¬ 
ting the initial value of the hydrogen absorption column 
density to the Galactic one of Nh = 3.4 x 10^'^ cm“^. We 
used the ’Neldermead’ optimisation method to And a global 
solution that resulted in a reduced x^/DoF of 1.14, a best-fit 
value of the photon index of F = 1.88±0.13, and a hydrogen 
absorption column density of Nh = (5±2) x 10^° cm“^. The 
resulting pile-up fraction was ~6 %, in agreement with our 
estimate above. 

To determine the flux we fit the spectra without a pile- 
up model and froze the power-law index to F = 1.88. We 
found an unabsorbed flux of (6.5±0.6) x 10”^^ erg cm~^ s~^ 
in the 0.3-10 keV range, corresponding to Lx = (9.0 ± 
0.8) X 10®® erg s“^ at the source distance of 3.39 Mpc. We 
also attempted to fit the same X-ray models to the Chandra 
spectrum as in Sec. 12.1.11 The spectrum is well fit by a 
power-law and does not require any additional component. 
The power law fit parameters and the X-ray luminosity are 
consistent with being the same for the Swift and the Chandra 
observations, which were sep a rated by ~4 months. 

Following IWalton et al.l ll2015ll , we fit the Chandra 
spectrum with two disc blackbody models in XSPEC 
(piLEUP*TBABS*TBABS*(DISKBB-fDISKBB)). We also fit a 
disc blackbody plus a p-free disk blackbody model 
(pileup*tbabs*tbabs*(diskbb-|-diskpbb)). These fits re¬ 
sulted in a reduced x^/DoF= 1.10 and x^/DoF= 1.09, re¬ 
spectively. These reduced x^ values may indicate that em¬ 
ploying two disks instead of a power law does not signifi¬ 
cantly improve the fit. The first hydrogen absorption column 
densities were frozen to the Galactic value. In both of the 
cases, the fit did not require any additional Nh on top of 
the Galactic one. The first fit resulted in inner disc temper¬ 
atures of Tin.i = 0.31 ± 0.04 keV and Tin ^2 = l-64o.2 keV. 
The second fit resulted in Tin,i = 0.28 ± 0.07 keV, Tin ,2 = 
2.1!to 7 keV, and p = O.OIq ®. The corresponding fitted tem¬ 
peratures between the first and the second fit agree within 
the errors. Furthermore, given that the value of p can be 
consistent with 0.75, with the standard disc, we cannot dif¬ 
ferentiate between a standard or a p-free disc model. 

The fitted parameters of the first fit are consistent with 
being the same for the Swift and the Chandra observations. 
However, the cooler compon ent is somewhat hotter than 
found by I Walton et al.l (l2015l l. On the other hand, if a disk 
blackbody plus p-free disk is employed, both, the cooler and 
hotter component a grees within the errors with the result of 
IWalton et al.l (120151') . 

2.2.2 VLA observations 

The Karl G. Jansky Very Large Array (VLA) observed 
Ho H X-1 simultaneously with Chandra on 25 May 2014. The 
observations were carried out under project code SF0330 
and ran from 03:44:50 UT for a total of 3.5 h. The VLA was 
in its most-extended A-configuration, and the correlator in¬ 
tegration time was 2 s. Two 1-GHz wide basebands covered 
the 8-10 GHz range (X band) with a total of 16 spectral 

^ http://cxc.harvard.edu/sherpa/threads/pileup/index.html 
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Table 2. Summary of results on the central component 


Instrument 

Band 

Flux density 

e 

Tb 



[hJy] 

[mas] 

[K] 

EVN 

1.6 GHz 

88±11 

16.4T1.7 

1.5x10® 

EVN 

5 GHz 

<36 

- 

- 

VLA 

8-10 GHz 

12±5 

<235 

- 

Instrument 

Band 

Unabsorbed flux 

r 

Nh 



[erg cm“^ 


[cm-2] 

Swift 

0.3-10 keV 

(6.6±0.4)x 10-12 

1.91T0.08 

(5±2)x1020 


Chandra 0.3-10 keV ( 6 . 5 ± 0 . 6 )XlO'^^ l.88±0.13 (5±2)xl02° 

The results of the Gaussian fits to the central radio component and the fit of an absorbed power-law model to the X-ray source. For the 

1.6-GHz VLBI data, the brightness temperature (Tu) value is also provided. 

be constrained at these signal-to-noise ratios. We modelled 
these components with elliptical Gaussians using the JM- 
FIT task within AIPS, and the results are shown in Ta¬ 
ble [S] We estimate the spectral index of the outer compo¬ 
nents between 5.24 and 9 GHz to be asw = —0.5 ± 0.4 and 
One = —0.8 ± 0.4, where Sv oc v°‘. However, we caution 
that the deconvolved size of the SW component is larger 
at X band than at C band. This could artificially flatten 
asw and the size of the components at X band are poorly 
constrained (Table |3]). 

For the central component, the size upper limits are 
in good agreement between the C-band uniformly-weighted 
image (< 0.28” x 0.19”) and the X-band robust 0 image (Ta¬ 
ble |3]). Assuming a non-variable source, the spectral index 
of the central component is ac = —4.7 ± 0.5. We discuss 
possible interpretations for such a steep spectrum in Sec. 
[33] Note that the minimum baseline length in the X-band 
data was 10.5 kA. To test whether resolution effects could be 
artificially steepening the spectrum, we produced a G-band 
image with a minimum baseline length of 10.5kA and found 
no change in either the peak brightness or the integrated flux 
density. This indicates that the unresolved central compo¬ 
nent at C band is unlikely to be contaminated by additional 
extended emission that was resolved out in our X-band ob¬ 
servations. 


o 

o 


70°42'16 



08 19 29 .6 
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Figure 2. The VLA image of Ho II X-1 in X band, using Briggs 
robust=l weighting. The lowest contours start at three times the 
rms noise level, at ±6.6 /rjy beam“^. Gontours are increased by 
factors of \/2. The peak brightness is 15 /ijy beam“^. The Gaus¬ 
sian restoring beam size is 0.41” X 0.22” at a major axis position 
angle of 75° E of N. The background colour image is arbitrarily 
scaled. 


windows, each with 64 2-MHz wide channels. We calibrated 
the data using CASA llMcMullin et al.|l2007f l version 4.1.0 
and the VLA pipeline version 1.2.0. The absolute amplitude 
and bandpass calibrator was 3C286 and the gain and phase 
calibrator was J0721±7120. The primary calibrator was ob¬ 
served once for 6 min and the secondary calibrator was ob¬ 
served for 50 s in every 7.5-min block. The total on-source 
time was 171 min. There was no self-calibration applied. Im¬ 
ages were made using Briggs robust=l and robust=0 weight¬ 
ing and the multi-frequency synthesis (MFS) method with 
two Taylor terms. The rms noise was 2.2 /ijy/beam and the 
synthesised beam FWHM was 0.41// x 0.22” in PA 75° (Fig. 

13). 

We detected three components, correspondi ng to those 
seen in C-band (4.5-6.5 GHz) VLA images llCseh et al.l 
[20 i 3). The central, NE, and SW components were detected 
with a signal to noise ratio of 6.6, 6.0, and 6.8, respectively. 
The in-band spectral index of the components could not 


3 DISCUSSION 


3.1 Is Ho II X-1 in the canonical X-ray hard 
state? 

It is common to fit ULX spectra with a two-component phe¬ 
nomenological model such as a power law and disk black- 
body. On the other hand, high-quality XMM-Newton spec¬ 
tra of ULXs sh ow a soft excess and a luminosity-dependen t 
spectral break llStobbart et al.l[2006l : iGladstoneetjaL 200^ 
and this is also the case for Ho H X-1 ( Kaiava et al.l 
l2012ll . that may indicate super-Eddington accret ion. On the 
other hand, ICaballero-Garcfa fc FabiarJ (l2010fl suggested 
that such a spectrum might be reflection-dominated and 
that the dominant contribution to the total luminosity 
is coming from the power-law component, potentially im¬ 
plying the sources to be sub-Eddington. However, in the 
case of Ho H X-1, this model resulted in a fitted pho¬ 
ton index that might be too steep for a canonical hard 




























6 D. Cseh et al. 


Table 3. The VLA X-band components 


Comp. 

Peak 

[/^Jy/b] 

Flux dens. 

[hJy] 

Size 

[mas X mas] 

Rob¬ 

ust 

PA 

Central 

12 ±3 

12 ±5 

< 230x240 

0 

76° 

Central 

13 ±2 

50± 10 

530±170x480±190 

1 

62° 

sw 

14 ±2 

64± 11 

700±140x450±110 

1 

91° 

NE 

12 ±2 

35 ±7 

500±160x300±200 

1 

94° 


The Gaussian fit results of the different components. 


state lICaballero-Garcfa fc FabiarJl201oh . Furthermore, NuS- 
TAR observations have confirmed th e presence of a high- 
energy cutoff for a number of ULXs llBachetti et al.l [20131 : 
IWalton et al .112013 : 1 Rana et al.ll2015h . which in those cases 
likely rules out a Compton hump arising from a reflection 
component in a canonical hard X-ray st ate. This has recent ly 
been shown specifically for Ho II X-1 dWalton et aLllioiSl l. 

Thus, despite the fact that the current spectra are well 
fit by a single power law, we argue that Ho H X-1 is unlikely 
to be in the canonical hard state. Additionally, there is no 
hint of a self-absorbed compact jet in Ho H X-1 (see later), as 
is typically seen in this state. It is also unlikel y that Ho H X- 
1 is in a non-standard, efficient hard state (iRushton et al.l 
bold: ICoriat et all l201lli since the X-rays are steady and 
decorrelated from the radio emission. Also, our X-ray spec¬ 
tra may as well be fit with two disk blackbodies. In this 
case, it would be straightforward to exclude Ho H X-1 being 
in a canonical hard state, and such a fit possibly indicates 
a non -standard high stat e or a super-Eddington accretion 
state ll Walton et al.l[2015l ~). 

In the following we investiga t e the temporal behaviour 
of the X-ray spectra. iGrise et al.l (1201011 presented multiple 
S'rui/f/XRT observations of Ho H X-1, finding the source 
to remain in a soft spectral state for over 4 months in 
2009/2010, in contrast to the hard spectra observed in 2006, 
which persisted for up to a month. The hard spectra from 
2006 could also be fit with a disk blackbody and a power law 
model. This is contrary to our spectra where this combina¬ 
tion of models resulted in a 0 keV inner disk temperature. 
On the other hand, the power-law components are identical 
within uncertainties between these Swift data sets, so a vari¬ 
able soft component is likely to be present. Also, comparing 
our Swift and Chandra spectra, the power-law components 
are consistent with being the same even though they are sep¬ 
arated by 4 months. Based on our observations, we conclude 
that the (apparently) hard spectra are in general steady, 
as compared to the soft spectra, which, beyond a possibly 
variable soft comp onent, are associated with irregular flares 
llCrise et aklbOKl ). The potential link between these flares 
and the radio activity needs further investigation. 


GHz. This non-detection at 5 GHz is consistent with the 
extended morphology found at 1.6 GHz, and also with the 
optically thin synchrotron spectrum observed at lower res¬ 
olution dCseh et al.ll201-i j. Thus, we infer that the VLBI 
emission is not due to a partially self-absorbed core jet. To 
estimate whether the source could be Doppler-boosted, we 
calculate its brightness temperature using: 

(l)-—o-a..,(^) (2^)^(==:)“ a, 

After substituting Sv = 88 /rjy at 1.66 GHz, 0 — 16.4 mas, 
and 2 = 0.00079, we find a brightness temperature of 
Tb ~ 1.5 X 10® K. This is a relatively low value with re¬ 
spect to the equipartitio n value of relativ istic compact jets 
of Tb,eq — 5 X 10^° K (iReadheadI Il994h . Given that the 
Doppler-factor is 5 ~ Tb/2b,eq, the low Tb of the emission 
argues against Doppler-boosted emission and hence against 
relativistic beaming being responsible for strong radio vari¬ 
ability (see next section). 

3.3 Time-dependent radio behaviour 

Assuming a non-variable source, the spectral index of the 
central component between 5.24 and 9.0 GHz is ac = —4.7± 
0.5, which is signific antly steeper tha n the previous estimate 
of a = —0.8 ± 0.2 dCseh et al.ll20ljj . In fact, such a steep 
spectrum is unphysical unless the synchrotron spectrum has 
a cutoff at the observed frequency due to, e.g., synchrotron 
ageing. We investigate this option in the following section. 

Since the two VLA measurements were taken about 
1.5 yr apart, we therefore consider time variability. Taking 
the 5.24-GHz VLA peak brightness of 152 /rjy/beam and a 
spectral index of a = —0.8 ± 0.2, the expected peak bright¬ 
ness of a point source at 9 GHz is 99±11 /rjy/beam. The 
detected X-band peak brightness of -^12 ^Jy/beam is at 
least a factor of 7.3 lower than this (a difference of ~ 8a). 
Given that resolution effects were excluded (Sec. 12. 2. 211 and 
that the central source is point-like at both C and X band, 
this indicates that the flux density decreased with time. 


3.2 VLBI morphology 

In light of the relatively stable X-ray spectra, in the follow¬ 
ing we shall investigate the contemporaneous behaviour of 
the central radio source that is coincident with the ULX 
position. 

The 1.6-GHz VLBI component was found to be ex¬ 
tended, with a size of ~0.26 pc. However, no milliarcsecond- 
scale radio emission was detected three months later at 5 


3.4 The jet properties of Ho II X-1 

The central component is due to optically thin synchrotron 
emission, consistent with a compact ejection. Since it was 
found to be a factor of ~3 brighter than the outer compo¬ 
nents, th is was interpreted as a renewed radio activity in 
the core dGseh et al.ll2014| j. The current VLA observation, 
obtained 1.5 yr later, also shows the same morphology. How¬ 
ever, the brightness of the central component is comparable 
to or less than that of the outer components. Therefore, the 
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central component has cooled much faster than the outer 
components. This is consistent with the deduced spectral 
indices in Sec. 12.2.21 The central component has apparently 
faded by at least a factor of 7.3, is resolved on VLBI scales, 
and is not affected by Doppler boosting. 

The variation in the flux density over the course of ~1.5 
yr is likely a result of expansion. On the other hand, with¬ 
out a light curve, and with a VLBI component that is not 
(yet) resolved into a two-sided jet, the evolutionary stage 
of the ejection is unclear. Here, we assume a spreading jet 
and a Mach num ber similar to that of the large-scale jet 
dCseh et al.ll20ljj . Then a lateral expansion velocity of 0.2c 
is expected for a speed of c along the jet axis. This max¬ 
imum expansion speed would result in a maximum radial 
increase of ~4 mas (8 mas FWHM) over a year and may 
limit the age of the ejection to ^2.1 years. On the other 
hand, the bulk speed is unknown and based on the geome¬ 
try it could be as slow as 0.17c. Considering such a mildly 
relativistic bulk speed, the radial expansion rate would be 
~1.3 mas/yr, implying an increased ejection age of -^13- 
46 yr for our measured size of 16.4-60 mas. 

In the following, we deduce an approximate adiabatic 
timescale (tad) and compare it with the synchrotron cool¬ 
ing timescale (r). By assuming a fixed expansion speed of 
Uexp = 0.2c, the adiabatic timescale can be approximated 
from the apparent physical size (r) of the ejection using 
tad — i—For an inner structure of size 16.4-60 mas, 

^ 'f^exp 

this gives tad — 3.2-11.7 yr. Next, we estimate the syn¬ 
chrotron lifetime of relativistic electrons at an observed fre¬ 
quency of 9 GHz, using r ~ 2.693 x (e.g. 

lLongaiilll994f) . where the units are in Hz, mG, and yr. We 
estimate the magnetic held strength from minimum energy 
arguments, assuming equipartition between the energy in 
relativistic particles and that in the magnetic held. Using 
B = 1.8 X VA mG le.g. lLongairj|l993 ~l. and 

substituting in the VLBI parameters of iz = 1.6 x 10® Hz, 
U = 3 X 10'*'^ m®, and U = 1.21 x 10^’’ WHz"\ we hnd 
B ~ 0 . 7777 ®^^ mG, where 77 — 1 is the ratio of energy in 
protons to that in relativistic electrons. Here we consider 
values of 77 = 1-100, and hnd r ~ (1.8-13) x 10® yi[f|. A 
comparison of these timescales indicates that the ejection 
is dominated by adiabatic losses, which results in the ob¬ 
served rapid hux decay. The outer components are bright 
because they are likely the working surfaces where the jet 
ejecta interact with the surroundings, allowing the conver¬ 
sion of kinetic energy to radiation. Finally, we note that 
with an average jet power of Qj ~ 2.1 x 10®® erg s“®, the 
surround ing radio bubble might have been inhated in just 
~390 yr JCseh et al.ll2014l L 


3.5 Summary and implications 

We have conducted quasi-contemporaneous X-ray and ra¬ 
dio observations of Ho H X-1. The X-ray observations show 
steady, apparently hard spectra with properties that are 
identical to the power-law component observed in spectra 


® Here we do not perform the calculations with the size upper 
limit of ~60 mas, because a higher size would mean even longer 
timescales. 


from 2006 llGrise et al.l 1201^ 1. By contrast, other observa¬ 
tions in 2009/2 010 have shown soft spectra and irregular 
flaring activity llGrise et al.ll2O10h . Given the steady, appar¬ 
ently hard emission, the X-ray flux variability may well arise 
purely from the soft component. The VLBI data reveal an 
extended radio source with a size of 0.26 pc. Based on its 
variability and brightness temperature, we find that it is 
consistent with non-thermal emission and that it is not af¬ 
fected by Doppler boosting. This is consistent with a sce¬ 
nario where the radio emission is due to a jet ejection event. 
Assuming a laterally-expanding jet, we estimate the dynam¬ 
ical age of the ejection to be ^2.1 yr. Whether an ejection 
event could be correlated with the soft or flaring nature of 
the X-ray spectra remains to be investigated. On the other 
hand, we argue that the hard X-ray spectra are unlikely 
to correspond to a canonical hard state. This conclusion is 
also supported by the lack of a classical flat-spectrum, par¬ 
tially self-absorbed compact radio jet. Also, the lack of any 
correlation between the X-rays and the radio emission may 
also rule out radiatively efficient hard states, i.e. the ‘outlier 
bran ch’ on the X-ray/ra dio correlation of BH X-ray binaries 
(e.g. iGoriat et al.l 1201111 . Using the VLA we have detected 
a triple structure at 8-10 GHz that corresponds well to the 
morphology seen at lower frequencies. However, the central 
component has faded by at least a factor of 7.3 in just 1.5 
years. We explain the observed cooling with adiabatic ex¬ 
pansion losses that dominate over synchrotron cooling, re¬ 
gardless of the jet composition (baryonic or leptonic). 

iGseh et al.l il2012l l reported unresolved radio emission in 
the vicinity of the ULX IG 342 X-1 and concluded that it 
was inconsistent with a self-absorbed compact jet, and could 
instead be a cl ump of emission from the radio nebula. On 
the other hand. [Marlowe et al.l (I2014! ) argued that this unre¬ 
solved emission could be consistent with an expanded ejec¬ 
tion similar to that in Ho H X-1. By comparison, it is indeed 
plausible that such a radio component in IG342 X-1 could 
have faded on a timescale of a few years. If t hat is the case, 
then the quoted size limit on the ejection dMarlowe et al.l 
l2014h of IG 342 X-1 might be too large, because the spa¬ 
tial variation would be small even with highly relativistic 
expansion speeds. 

In light of our results on Ho H X-1, it is possible that 
other ULX radio bubbles could also have been inflated by 
transient ejecta. This could imply that it is not necessary to 
invoke self-absorbed compact jets to inflate ULX bubbles. 
Moreover, the self-absorbed compact jet in Gyg X-1 has been 
argued to be incapable of inflating its bubble llHein 3 l2006ll 
due to the mismatch between the energetics of the jet and 
the bubble. On the other hand, expanding ejecta and disk 
winds, which are potentially triggered by high accretion-rate 
events, could serve as a more plausible explanation in the 
above cases, as they might indeed be capable of carrying the 
required energy. 

Ho H X-1 shows a collimated jet structure that is cou¬ 
pled with a high average X-ray luminosity whose ionizing 
rate is comparable to the average jet power. This may be 
unexpected given that AGNs with high acc retion rate s show 
strong winds but n o powerful jets (e.g. I Gan et ^ I I20i 4 
iTombesi et al.l[2015l L Similarly, in the context of regular BH 
X-ray bin aries the jet is also suppressed in disc-dominated 
states, and lMei^ (Il996l l argues that a super-Eddington wind 
could be a secondary jet-suppression mechanism. Nonethe- 
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less, Ho II X-1 does show discrete jet ejecta. This may imply 
that if the current X-ray luminosity is super-Eddington, then 
Ho H X-1 has to switch between sub- and super-Eddington 
rates to form jet ejecta. Alternatively, it could be that the 
source is not super-Eddington, but is rather in a persistent 
high, intermediate state where jet ejecta are thought to be 
released whenever crossing the so-called ‘jet line’. 

X-ray monitoring of the source with triggered radio ob¬ 
servations could differentiate between any of the cases and 
test if X-ray flaring traces a jet ejection event. 
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